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Abstract. The present experimental status in the search for neutrinoless double beta 
decay is reviewed, with emphasis on the first indication for neutrinoless double beta 
decay found in the HEIDELBERG-MOSCOW experiment, giving first evidence for 
lepton number violation and a Majorana nature of the neutrinos. Future perspectives 
of the field are briefly outlined. 



1 Introduction 

The neutrino oscillation interpretation of the atmospheric and solar neutrino 
data, deliver a strong indication for a non-vanishing neutrino mass. While such 
kind of experiments yields information on the difference of squared neutrino mass 
eigenvalues and on mixing angles (for the present status see, e.g. [50,51]), the 
absolute scale of the neutrino mass is still unknown. Information from double 
beta decay experiments is indispensable to solve these questions [5,6]. Another 
important problem is that of the fundamental character of the neutrino, whether 
it is a Dirac or a Majorana particle [7,8]. Neutrinoless double beta decay could 
answer also this question. Perhaps the main question, which can be investigated 
by double beta decay with high sensitivity, is that of lepton number conservation 
or non-conservation. 

Double beta decay, the rarest known nuclear decay process, can occur in 
different modes: 

2i//3/3 - decay : A(Z, N) -> A(Z+2, N-2) + 2e~ + 2P e (1) 

Ou/3/3 - decay : A{Z, N) -> A(Z+2, N-2) + 2e~ (2) 

0u(2) x P(3 - decay : A(Z, N) -> A(Z+2, N-2) + 2e~ + (2)x (3) 



While the two- neutrino mode (1) is allowed by the Standard Model of particle 
physics, the neutrinoless mode (OV/3/3 ) (2) requires violation of lepton number 
(AL=2). This mode is possible only, if the neutrino is a Majorana particle, i.e. 
the neutrino is its own antiparticle (E. Majorana [7], G. Racah [8], for subsequent 
works we refer to [10-12], for some reviews see [24-27,6,28]). First calculations 
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of Qvj3j3 decay based on the Majorana theory have been done by W.H. Furry [9]. 
The most general Lorentz-invariant parametrization of neutrinoless double beta 
decay is shown in Fig. 1. 




(c) (d) 



Fig. 1. Feynman graphs of the general double beta decay with long range (a-c) and 
short range (d) parts (see [14]). 

The usually used assumption is that the first term (i.e. the Majorana mass 
mechanism) dominates the decay process. However, as can be seen from Fig. 1, 
and as discussed elsewhere (see, e.g. [6,29,30,42]) neutrinoless double beta decay 
can not only probe a Majorana neutrino mass, but various new physics scenarios 
beyond the Standard Model, such as R-parity violating supersymmetric models, 
R-parity conserving SUSY models, leptoquarks, violation of Lorentz-invariance, 
and compositeness (for a review see [6,29,30]). Any theory containing lcpton 
number violating interactions can in principle lead to this process allowing to 
obtain information on the specific underlying theory. It has been pointed out 
already in 1982, however, that independently of the mechanism of neutrinoless 
double decay, the occurence of the latter implies a non-zero neutrino mass and 
vice versa [13]. This theorem has been extended to supersymmetry. It has been 
shown [15] that if the neutrino has a finite Majorana mass, then the sncutrino 
necessarily has a (B-L) violating 'Majorana' mass, too, independent of the mech- 
anism of mass generation. The experimental signature of the neutrinoless mode 
is a peak at the Q-value of the decay. 

Restricting to the Majorana mass mechanism, a measured half-life allows to 
deduce information on the effective Majorana neutrino mass (to) , which is a 
superposition of neutrino mass eigenstates [24,25]: 

[^°; 2 (0+ - 0+)]- 1 = C mm ^L- + C m { v ) 2 + C\\ (A) 2 + C mn (r,)^- 

+C mX (X)^+C vX (r 1 )(X), (1) 
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( m ) = | m (l)| +e *| m (2)| + e ^3| m (3)| ( (2) 

where mil = \fnee | exp (ifa) (i = 1,2,3) arc the contributions to (m) from in- 
dividual mass eigenstates, with fa denoting relative Majorana phases connected 
with CP violation, and C m m,C mi ... denote nuclear matrix elements squared, 
which can be calculated, (see, e.g. [33], for a review and some recent discussions 
see e.g. [34,6,25,26,37,40-42]). Ignoring contributions from right-handed weak 
currents on the right-hand side of eq.(l), only the first term remains. 

The effective mass is closely related to the parameters of neutrino oscillation 
experiments, as can be seen from the following expressions 

\m£\ = \U el \ 2 m u (3) 
\m^\ = \U e2 \ 2 ^Am 2 1+ m 2 , (4) 
\mi 3 J\ = \U e3 \ 2 ^Am 2 2 +Am 2 21 +ml (5) 

Here, U e i are entries of the neutrino mixing matrix, and Amfj = |m| — m 2 \, 
with rm denoting neutrino mass eigenstates. U e i and Am 2 can be determined 
from neutrino oscillation experiments. 

The importance of (m) for solving the problem of the structure of the neutrino 
mixing matrix and in particular to fix the absolute scale of the neutrino mass 
spectrum which cannot be fixed by v - oscillation experiments alone, has been 
discussed in detail in e.g. [3,5,43]. 

Double beta experiments to date gave only upper limits for the effective 
mass. The most sensitive limits [17-19] were already of striking importance for 
neutrino physics, excluding for example, in hot dark matter models, the small 
mixing angle (SMA) MSW solution of the solar neutrino problem [44-49,5,6] in 
degenerate neutrino mass scenarios. 

The HEIDELBERG-MOSCOW double beta decay experiment in the Gran 
Sasso Underground Laboratory [2,21,23,29,22,6] searches for double beta decay 
of 76 Ge — > 76 Se + 2 e~ + (2D) since 1990. It is the most sensitive double beta 
experiment since almost eight years now. The experiment operates five enriched 
(to 86%) high-purity 76 Ge detectors, with a total mass of 11.5 kg, the active 
mass of 10.96 kg being equivalent to a source strength of 125.5 mol 76 Ge nuclei. 

In this paper, we present a new, refined analysis of the data obtained in the 
HEIDELBERG-MOSCOW experiment during the period August 1990 - May 
2000 which have recently been published [19]. The analysis concentrates on the 
neutrinoless decay mode which is the one relevant for particle physics (see, e.g. 
[6]). First evidence for the neutrinoless decay mode will be presented (a short 
communication has been given already in [1] (see also [105]), and first reactions 
have been published already [3,52-75]). For the results concerning 2v(3(3 decay 
and Majoron-accompanied decay we refer to [19]. The results will be put into 
the context of other present (3(3 activities. 
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2 Experimental Set-up and Results 

A detailed description of the HEIDELBERG-MOSCOW experiment has been 
given recently in [16], and in [19]. Therefore only some important features will 
be given here. We start with some general notes. 

1. Since the sensitivity for the Qv(3(3 half- life is 

(and /*ov ~ ( m ^))j w hh a denoting the degree of enrichment, e the efficiency of 
the detector for detection of a double beta event, M the detector (source) mass, 
AE the energy resolution, B the background and t the measuring time, the sen- 
sitivity of our 11 kg of enriched 76 Ge experiment corresponds to that of an at last 
1.2 ton natural Gc experiment. After enrichment, energy resolution, background 
and source strength are the most important features of a (3(3 experiment. 

2. The high energy resolution of the Ge detectors of 0.2% or better, assures 
that there is no background for a 0v(3(3 line from the two-neutrino double beta 
decay in this experiment, in contrast to most other present experimental ap- 
proaches, where limited energy resolution is a severe drawback. 

3. The efficiency of Ge detectors for detection of Qv(3(3 decay events is close 
to 100% (95%, see [76]). 

4. The source strength in this experiment of 11kg is the largest source 
strength ever operated in a double beta decay experiment. 

5. The background reached in this experiment is, with 0. 17 events/kg ykeV 
in the Qv(3(3 decay region (around 2000-2080 keV), the lowest limit ever obtained 
in such type of experiment. 

6. The statistics collected in this experiment during 10 years of stable running 
is the largest ever collected in a double beta decay experiment. 

7. The Q value for neutrinoless double beta decay has been determined re- 
cently with very high precision to be 0^=2039.006(50) keV [77,78]. 

We give now some experimental details. All detectors (whose technical pa- 
rameters are given in Table 1 (see [16])), except detector No. 4, are operated 





Total 


Active 


Enrichment 


FWHM 1996 


FWHM 2000 


Detector 


Mass 


Mass 


in '"Ge 


at 1332 keV 


at 1332 keV 


Number 


[kg] 


[kg] 


[%] 


[keV] 


[keV] 


No. 1 


0.980 


0.920 


85.9 ± 1.3 


2.22 ± 0.02 


2.42 ± 0.22 


No. 2 


2.906 


2.758 


86.6 ± 2.5 


2.43 ± 0.03 


3.10 ± 0.14 


No. 3 


2.446 


2.324 


88.3 ± 2.6 


2.71 ± 0.03 


2.51 ± 0.16 


No. 4 


2.400 


2.295 


86.3 ± 1.3 


2.14 ± 0.04 


3.49 ± 0.24 


No. 5 


2.781 


2.666 


85.6 ± 1.3 


2.55 ± 0.05 


2.45 ± 0.11 



Table 1. Technical parameters of the five enriched detectors. (FWHM - fall width at 
half maximum.) 
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Fig. 2. The /3/?-laboratory of the HEIDELBERG-MOSCOW experiment in the Gran 
Sasso (upper left) and its location between halls A and B (upper right), and four of 
the enriched detectors during installation (lower part). 

in a common Pb shielding of 30 cm, which consists of an inner shielding of 10 
cm radiopure LC2-grade Pb followed by 20 cm of Boliden Pb. The whole setup 
is placed in an air-tight steel box and flushed with radiopure nitrogen in order 
to suppress the 222 Rn contamination of the air. The steel box is centered inside 
a 10-cm boron-loaded polyethylene shielding to decrease the neutron flux from 
outside. An active anticoincidence shielding is placed on the top of the setup to 
reduce the effect of muons. Detector No. 4 is installed in a separate setup, which 
has an inner shielding of 27.5 cm electrolytical Cu, 20 cm lead, and boron-loaded 
polyethylene shielding below the steel box, but no muon shielding. Fig. 2 gives 
a view of the experimental setup. 

To check the stability of the experiment, a calibration with a 228 Th and a 
152 Eu+ 22S Th, and a 60 Co source is done weekly. High voltage of the detectors, 
temperature in the detector cave and the computer room, the nitrogen flow 
in the detector boxes, the muon anticoincidence signal, leakage current of the 
detectors, overall and individual trigger rates are monitored daily. The energy 
spectrum is taken in 8192 channels in the range from threshold up to about 3 
MeV, and in a parallel spectrum up to about 8 MeV. 

Because of the big peak-to-Compton ratio of the large detectors, external 7 
activities are relatively easily identified, since their Compton continuum is to a 
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large extent shifted into the peaks. The background identified by the measured 
7 lines in the background spectrum consists of: 

1. primordial activities of the natural decay chains from 238 U, 232 Th, and 40 K 

2. anthropogenic radio nuclides, like 137 Cs, 134 Cs, 125 Sb, 207 Bi 

3. cosmogenic isotopes, produced by activation due to cosmic rays. 

The activity of these sources in the setup is measured directly and can be 
located due to the measured and simulated relative peak intensities of these 
nuclei. 

Hidden in the continuous background are the contributions of 

4. the bremsstrahlungs spectrum of 210 Bi (daughter of 210 Pb), 

5. elastic and inelastic neutron scattering, and 

6. direct muon-induced events. 

External a and (3 activities are shielded by the 0.7- mm inactive zone of the p-type 
Gc detectors on the outer layer of the crystal. The enormous radiopurity of HP- 
germanium is proven by the fact that the detectors No. 1, No. 2 and No. 3 show 
no indication of any a peaks in the measured data. Therefore no contribution of 
the natural decay chains can be located inside the crystals. Detectors No. 4 and 
No. 5 seem to be slightly contaminated with 210 Pb on the level of few /xBq/kg, 
most likely surface contaminations at the inner contact. This contamination was 
identified by a measured a peak in the background spectrum at 5.305 MeV of the 
daughter 210 Po and the constant time development of the peak counting rate. 
There is no contribution to the background in the interesting evaluation areas 
of the experiment due to this activity. For further details about the experiment 
and background we refer to [16,76] (see also Table 2). 

In the vicinity of the Q- value of the double beta decay of Qpp = 2039.006(50) 
keV, very weak lines at 2034.744 and 2042 keV from the cosmogenic nuclide 56 Co, 
and from 21A Bi ( 238 J7-decay chain) at 2010.7, 2016.7, 2021.8 and 2052.9 keV, may 
be expected. 

On the other hand, there are no background 7-lines at the position of an ex- 
pected OvPfl line, according to our Monte Carlo analysis of radioactive impurities 
in the experimental setup [76] and according to the compilations in [79] . 

In total 55 possible 7-lines from various isotopes in the region between 2037 
and 2041 keV are known [79]. Only 5 of the isotopes responsible for them ( 102 Rh, 
146 Eu, 124 I, 124 Sb and 170 Lu) have half-lifes larger than 1 day. However, some of 
the isotopes yielding lines in this energy range can in principle be produced by 
inelastic hadron reactions (induced by muons or neutrons). 

Therefore each of these 55 isotopes was checked for the existence of a 7-line 
from the isotope which has a high emission probability Ih ■ A search was made for 
this 7-line in the measured spectrum to obtain its intensity (Sh) or an upper limit 
for it. Then the adopted intensity So for a 7-line from the same isotope in the 
area around ~2039keV can be calculated by using the emission probability I 
for the line at ^2039 keV. Different absorption for gammas of different energies 
are taken into account in a schematic way. 

If the calculated intensity Si < 1 for the 7-line in the interesting area, this 
isotope can be safely excluded to contribute a significant part to the background. 
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Detec- 
tor 


Life 
Time 


Date 


Shielding 


Background*) 
[counts/ 
keV y kg] 


PSA 


Number 


[ days ] 


Start End 


Cu 


Pb 


boron- 
poly. 


2000.-2100. 

keV 




No.l 


387.6 


8/90-8/91 


yes 


0.56 


no 






1/92-8/92 






no 


No.2 


225.4 


9/91 - 8/92 


yes 


0.29 


no 






Common 


shielding lor three detectors 




No.l 


382.8 


9/92 - 1/94 


yes 


0.22 


no 


No.2 


383.8 


9/92 - 1/94 


yes 


0.22 


no 


No.3 


382.8 


9/92 - 1/94 


yes 


0.21 


no 


No.l 


263.0 


2/94 - 11/94 




yes 


yes 


0.20 


no 


No.2 


257.2 


2/94 - 11/94 




yes 


yes 


0.14 


no 


No.3 


263.0 


2/94 - 11/94 




yes 


yes 


0.18 


no 


Full Setup 

Four detectors in common shielding, one detector separate 


No.l 


203.6 


12/94 - 8/95 




yes 


yes 


0.14 


no 


No.2 


203.6 


12/94 - 8/95 




yes 


yes 


0.17 


no 


No.3 


188.9 


12/94 - 8/95 




yes 


yes 


0.20 


no 


No.5 


48.0 


12/94 - 8/95 




yes 


yes 


0.23 


since 2/95 


No.4 


147.6 


1/95 - 8/95 


yes 


0.43 


no 


No.l 


203.6 


11/95 - 05/00 




yes 


yes 


0.170 


no 


No.2 


203.6 


11/95 - 05/00 




yes 


yes 


0.122 


yes 


No.3 


188.9 


11/95 - 05/00 




yes 


yes 


0.152 


yes 


No.5 


48.0 


11/95 - 05/00 




yes 


yes 


0.159 


yes 


No.4 


147.6 


11/95 - 05/00 


yes 


0.188 


yes 



Table 2. Development of the experimental set-up and of the background numbers in 
the different data acquisition periods for the enriched detectors of the HEIDELBERG- 
MOSCOW experiment. *) Without PSA method. 



For example, the isotope 39 Xe possesses a 7-line at 2039.1 keV with an emis- 
sion probability of 0.078%. This isotope also possesses a 7-line at 225.4 keV with 
an emission probability of 3.2%. The intensity of the line at 225.4 keV in our spec- 
trum was measured to be <6.5 counts. This means that the 7-line at 2039.1 keV 
has < 2 ° 3 °^ 8 6.5 = 0.32 counts, and therefore can be excluded. 

Only eight isotopes could contribute a few counts according to the calculated 
limits, in the interesting area for the Qvj3(3 -decay area: 52m Fe, 93m Ru, 120 In, 
131 Ce, 170 Lu, 170m Ho, m Ta and 198 T1. Most of them have a half-life of a few 
seconds, only 170 Lu has a half-life of 2.01 days, and no one of them has a longer 
living mother isotope. To contribute to the background they must be produced 
with a constant rate, e.g. by inelastic neutron and/or muon reactions. Only 5 
isotopes can be produced in a reasonable way, by the reactions listed in Tabl. 3. 

Except 120 Sb each of the target nuclides is stable. All reactions induced 
with a-particles can be excluded due to the very short interaction length of 
a-particles. Two possibilities remaining to explain possible events in the 0^/3/3 
-decay area would be: 
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isotope 


52m Fe 


M;to Ru 


i2u In 


17Um Ho 


l«8 T j 


production 
reaction 


50 Cr(Q,2n 7 ) 


92 Mo(a,3n) 


120 Sn(n,p), 
123 Sb(n,a) 


17 °Er(n,p) 


197 Au(q, 3n 7 ) 



Table 3. Possible reactions (second line) which could produce the isotopes listed in 
the first line (see [83]). 



. 120 Sn(n,p) 120 In: 

The cross section for this reaction is 2.5±1 mb for E n = 14.5 MeV [84]. 
Assuming a neutron flux of (0.4±0.4) x 10~ 9 em \s s for neutrons with an energy 
between 10-15 MeV as measured in the Gran Sasso [85] the rate of 120 In 
atoms produced per year is about 2 x 10~ 5 when there are 50 g of 120 Sn in 
the detector setup. Even when the cross-section is larger for lower energies, 
this can not contribute a significant number of counts to the background. 
• 170 Er(n,p) 170m Ho: 

The cross section for this reaction is about 1.13±lmb for E n — 14.8 MeV 
[86]. Assuming again a neutron flux of (0.4±0.4)xl0~ 9 em \; - for neutrons 
with an energy between 10-15 MeV the rate of 170 Er atoms produced per 
year is even less when assuming 50 g of 120 Er in the detector-setup. 

In both cases it would be not understandable, how such large amounts of 
120 Sn or 170 Er could have come into the experimental setup. Concluding we do 
not find indications for any nuclides, that might produce 7-lines with an energy 
around 2039 keV in the experimental setup. 

Fig. 3 shows the combined spectrum of the five enriched detectors obtained 
over the period August 1990 - May 2000, with a statistical significance of 54.981 
kgy (723.44 molyears). (Note that in Fig. 1 of [19]) only the spectrum of the first 
detector is shown, but normalized to 47.4 kgy [20]). The identified background 
lines give an indication of the stability of the electronics over a decade of mea- 
surements. The average rate (sum of all detectors) observed over the measuring 
time, has proven to be constant within statistical variations. Fig. 4 shows the 
part of the spectrum shown in Fig. 3, in more detail around the Q-value of double 
beta decay. Fig. 5 shows the spectrum of single site events (SSE) obtained for 
detectors 2,3,5 in the period November 1995 - May 2000, under the restriction 
that the signal simultaneously fulfills the criteria of all three methods of pulse 
shape analysis we have recently developed [80,82] and with which we operate all 
detectors except detector 1 (significance 28.053 kgy) since 1995. 

Double beta events are single site events confined to a few mm region in the 
detector corresponding to the track length of the emitted electrons. In all three 
methods, mentioned above, the output of the charge-sensitive preamplifiers was 
differentiated with 10-20 ns sampled with 250 MHz and analysed off-line. The 
methods differ in the analysis of the measured pulse shapes. The first one relies 
on the broadness of the charge pulse maximum, the second and third one are 
based on neural networks. All three methods are 'calibrated' with known double 
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energy [keV] 



Fig. 3. Sum spectrum of enriched detectors Nr. 1,2,3,4,5 over the period August 1990 - 
May 2000 (54.9813 kgy) measured in the HEIDELBERG-MOSCOW experiment (bin- 
ning 0.36 keV). The sources of the main identified background lines are noted. 

escape (mainly SSE) and total absorption (mainly MSE) 7-lines [80-82,76]. They 
allow to achieve about 80% detection efficiency for both interaction types. 

The expectation for a 0v(3j3 signal would be a line of single site events on 
some background of multiple site events but also single site events, the latter 
coming to a large extent from the continuum of the 2614 keV 7-line from 208 TZ 
(see, e.g., the simulation in [80]). From simulation we expect that about 5% of 
the double beta single site events should be seen as MSE. This is caused by 
bremsstrahlung of the emitted electrons [76] . 

Installation of PSA has been performed in 1995 for the four large detectors. 
Detector Nr. 5 runs since February 1995, detectors 2,3,4 since November 1995 
with PSA. The measuring time with PSA from November 1995 until May 2000 
is 36.532kg years, for detectors 2,3,5 it is 28.053 kgy. 

Fig. 6 shows typical SSE and MSE events from our spectrum. 
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2060 2070 2080 

energy [keV] 

Fig. 4. Sum spectrum of the 76 Ge detectors Nr. 1,2,3,4,5 over the period August 1990 
to May 2000, (54.9813 kgy) in the energy interval 2000 - 2080 keV, around the Q w 
value of double beta decay (Qpp = 2039.006(50) keV) summed to IkeV bins. The 
curve results from Bayesian inference in the way explained in sec. 3. It corresponds to 
a half-life T^ 2 = (0.80 - 35.07) x 10 25 y (95% c.l.) 




2060 2070 2080 

energy [keV] 

Fig. 5. Sum spectrum of single site events, measured with the detectors Nr. 2,3,5 oper- 
ated with pulse shape analysis in the period November 1995 to May 2000 (28.053 kgy), 
summed to 1 keV bins. Only events identified as single site events (SSE) by all three 
pulse shape analysis methods [80-82] have been accepted. The curve results from 
Bayesian inference in the way explained in sec. 3. When corrected for the efficiency 
of SSE identification (see text), this leads to the following value for the half-life: 
T°^ 2 =(0.88 - 22.38)x 10 25 y (90% c.l.). 



All the spectra are obtained after rejecting coincidence events between dif- 
ferent Ge detectors and events coincident with activation of the muon shield. 
The spectra, which are taken in bins of 0.36 keV, are shown in Figs. 4,5, Fig. 2 of 
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Fig. 6. Left: Shape of one candidate for 0i/f3/3 decay (energy 2038.61 keV) classified as 
SSE by all three methods of pulse shape discrimination. Right: Shape of one candidate 
(energy 2038.97 keV) classified as MSE by all three methods. 



[1] in 1 keV bins, which explains the broken number in the ordinate. We do the 
analysis of the measured spectra with (Fig. 4) and without the data of detec- 
tor 4 (see Fig. 2 in [1], 46.502 kg y) since the latter does not have a muon shield 
and has the weakest energy resolution. The 0.36 keV bin spectra are used in all 
analyses described in this work. We ignore for each detector the first 200 days of 
operation, corresponding to about three half-lives of 56 Co (7\/ 2 = 77.27 days), 
to allow for some decay of short-lived radioactive impurities. 



a 1000 



800 



600 



400 



200 




1430 



1440 



1460 



1470 



1480 1490 
Energy [keV] 



Fig. 7. Fit of one the lines (here the 1460.81 keV line from 40 K) in the sum spectrum of 
all five detectors, used for the calibration and the determination of the energy resolution 
(see Table 4). 



The background rate in the energy range 2000 - 2080 keV is found to be (0.17 
± 0.01)events/kgykeV (without pulse shape analysis) considering all data as 
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background. This is the lowest value ever obtained in such type of experiments. 
The energy resolution at the Qpp value in the sum spectra is extrapolated from 
the strong lines in the spectrum to be (4.00 ± 0.39) keV in the spectra with 
detector 4, and (3.74 ± 0.42) keV (FWHM) in the spectra without detector 
4 (see Fig. 7 and Table 4). The energy calibration of the experiment has an 
uncertainty of 0.20 keV (see Table 4). 



energy [keV] 


energy [keV] 


width [keV] 


width [keV] 


fit 


from [79] 


fit 


from calc. 


1460.81 ± 0.02 


1460.81 


1.49 ± 0.01 


1.49 ± 0.13 


1764.56 ± 0.05 


1764.49 


1.70 ± 0.05 


1.59 ± 0.15 


2103.31 ± 0.45 


2103.53 


1.86 ± 0.35 


1.71 ± 0.16 


2204.12 ± 0.14 


2204.19 


1.89 ± 0.13 


1.74 ± 0.17 


2447.73 ± 0.26 


2447.86 


1.82 ± 0.33 


1.82 ± 0.18 


2614.48 ± 0.07 


2614.53 


1.80 ± 0.06 


1.88 ± 0.18 



Table 4. Energies and widths of some prominent lines in the sum spectrum of all five 
detectors determined by our energy calibration and peak fit methods, and comparison 
with the energy given in the literature [79], and the fitted dependence of the width as 
function of energy. 



3 Data Analysis 

We analyse the measured spectra with the following methods: 

1. Bayesian inference, which is used widely at present in nuclear and as- 
trophysics (see, e.g. [87,89,91,108]). This method is particularly suited for low 
counting rates, where the data follow a Poisson distribution, that cannot be 
approximated by a Gaussian. 

2. Feldman-Cousins Method [88,89]. 

3. Maximum Likelihood Method (see [89,92]). 

The Bayesian method is described in the next subsection, 3.1. In subsection 
3.2 we give some numerical examples of the sensitivities of this method and of 
the Maximum Likelihood Method in the search for events of low statistics [108]. 

3.1 The Bayesian Method 

We first describe the procedure summarily and then give some mathematical 
details (see [108]). 

One knows that the lines in the spectrum are Gaussians with the standard 
deviation a=1.70keV in Fig. 4 and cr=1.59keV in Fig. 5. This corresponds to 
4.0(3.7) keV FWHM. Given the position of a line, we used Bayes theorem to 
infer the contents of the line and the level of a constant background. 

Bayesian inference yields the joint probability distribution for both param- 
eters. Since we are interested in the contents of the line, the other parameter 
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was integrated out. This yields the distribution of the line contents. If the distri- 
bution has its maximum at zero contents, only an upper limit for the contents 
can be given and the procedure does not suggest the existence of a line. If the 
distribution has its maximum at non-zero contents, the existence of a line is sug- 
gested and one can define the probability that there is a line with non-zero 
contents. 

We define the Bayesian procedure in some more detail. It starts from the 
distribution p(xi...xm\p,v) °f the count rates x\...Xm in the bins 1..M of the 
spectrum - given two parameters p, rj. The distribution p is the product 



M 



p(x!...x m \p,v) = IT — r e " 



(7) 



of Poissonians for the individual bins. The expectation value is the superpo- 
sition 

A fc = p[ V fi(k) + (l- V ) f 2 (k)] (8) 

of the form factors fi of the line and f 2 of the background; i.e. fi(k) is the 
Gaussian centered at E with the above-mentioned standard deviation value and 
f 2 (k)=f 2 is a constant. Note that the model allows for a spectrum of background 
only, i.e. i]—0, and in this sense also tests the hypothesis 'only background'. 
Since 



M 



for v =1,2, 



fc=i 



one has 



A I 



k=l 



(9) 



(10) 



Hence, p parametrizes the total intensity in the spectrum, and 77 is the relative 
intensity in the Gaussian line. 

The total intensity p shall be integrated out. For this purpose, one needs the 
prior distribution p{p\ri) of p for fixed r\. We obtain it from Jeffreys' rule (§5.35 
of [91]). 

1/2 

p(p\v) = 



Q2 



dp 2 



\np 



(11) 



The overline denotes the expectation value with respect to x\...xu- 
The integration 

p 1 (x 1 ...x M \v) ~ J dpp(x 1 ...x M \p,ri) p(p\v) 



(12) 



then yields the model pi conditioned by r\ alone. It is normalized to unity and 
the prior distribution 

_ 1/2 



Pi (v) 



Q2 

drj 



2 lnpi 



(13) 
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of r] is obtained by application of Jeffrey's rule to pi. 1 Bayes' theorem yields the 
posterior distribution 

□ i i \ pi (xi...x M \v) Mn) nA , 

P 1 {r]\x 1 ...XM) = ri ( 14 ) 

J drip! fix 

of 77. 

From the posterior the 'error interval' for 77 is obtained. It is the shortest 
interval in which 77 lies with probability K. The length of an interval is defined by 
help of the measure yi\(r\). We call this the Bayesian interval for the probability 
K in order to distinguish it from a confidence interval of classical statistics. There 
is a limit, where Bayesian intervals agree with confidence intervals. See below. 

The borders of a Bayesian interval are given by the intersections of the like- 
lihood function P\ (j]\x) j ' [i\{r\) ~ pi(a;|r7) with a horizontal line at rji, rjh (see 
Fig. 8). The probability K is obtained by integrating Pi from r\i to r\h- 



3 p-, — , — , — , — 1 — , — , — , — , — 1 — , — , — , — , — 1 — , — , — , — , — 1 — , — , — , — ^ 




Fig. 8. The figure shows the relation between Pi, K and Kb. K is the integral over Pi 
in the limits [f?i,»7h]. The integral over Pi in the interval [0,t]e] is Kg. 



When the likelihood function has its maximum at 77 = 0, then the Bayesian 
interval will - for every K - include 77 = 0. Then this value cannot be excluded 
and only an upper limit for the contents of the line can be given. 

When the maximum of the likelihood function is at a point different from 
77 = - as it is in Fig. 8- then there is a range of K- values such that the associated 
interval excludes the point 77 = 0. Under this condition let us construct the 
interval that has its lower border at 77 = 0. It extends up to 7]e ■ The associated 
probability is called Ke- The point 77 = now limits the possible 77- values in a 

1 We have done the analysis (sections 3.3, 3.4) also using other prior distributions and 
found only little effect on the result (see [106]). 
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non-trivial way because for every K < Kg , the associated error interval excludes 
zero. We call Ke the probability that there exists a line. 

The above considerations lead to a peak finding procedure [108]. One can 
prescribe a line at an arbitrary energy E of the spectrum - say of Fig. 4 - and 
determine the probability Kg that there is one. Such searches lead to the results 
given in the next section. 

Let us note that classical and Bayesian statistics become equal to each other 
when the likelihood function is well approximated by a Gaussian. In this case, 
the probability K is the same as classical confidence. 

Note that the method of minimum % 2 is based on an even more stringent 
limit. It requires Gaussian distributions of the data. Since the Gaussian is de- 
fined everywhere on the real axis, the method can yield negative values of the 
parameters, especially negative rj in the present case. The Bayesian method re- 
spects the natural limitations of the parameters because it accepts non-Gaussian 
distributions. 

The method of maximum likelihood is, roughly spoken, the Bayesian method 
with the prior distribution set constant. This is a useful approximation when 
the posterior is sufficiently narrow. Then the posterior becomes approximately 
Gaussian. In this sense, the method of maximum likelihood is based on a hidden 
Gaussian approximation. 

3.2 Numerical Simulations 

To check the methods of analysing the measured data (Bayes and Maximum 
Likelihood Method), and in particular to check the programs we wrote, we have 
generated spectra and lines with a random number generator and performed then 
a Bayes and Maximum Likelihood analysis (see also [107]). The length of each 
generated spectrum is 8200 channels, with a line located at bin 5666, the width of 
the line (sigma) being 4 channels (These special values have been choosen so that 
every spectrum is analogue to the measured data). The creation of a simulated 
spectrum is executed in two steps, first the background and second the line was 
created, using random number generators available at CERN (see [93]). In the 
first step, a Poisson random number generator was used to calculate a random 
number of counts for each channel, using a mean value of n=A or /i = 0.5 , 
respectively, in the Poisson distribution 

n 

P( n ) = ft- e-" (15) 

These mean values correspond roughly to our mean background measured in the 
spectra with or without pulse shape analysis. 

In the second step, a Gaussian random number generator was used to calcu- 
late a random channel number for a Gaussian distribution with a mean value of 
5666 (channel) and with a sigma of 4 (channels). The contents of this channel 
then is increased by one count. This Gaussian distribution filling procedure was 
repeated for n times, n being the number of counts in this line. 
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For each choice of n and n, 100 different spectra were created, and analysed 
subsequently with two different methods: the maximum likelihood method (using 
the program set of [92]) and the Bayes-method. Each method, when analysing a 
spectrum, gives a lower and an upper limit for the number of counts in the line 
for a given confidence level c (e.g. c = 95%) (let us call it confidence area A). 
A confidence level of 95% means, that in 95% of all cases the true value should 
be included in the calculated confidence area. This should be exactly correct 
when analysing an infinite number of created spectra. When using 100 spectra, 
as done here, it should be expected that this number is about the same. Now 
these 100 spectra with a special n and fi are taken to calculate a number d, which 
is the number of that cases, where the true value n is included in the resulting 
confidence area A. 

This number d is given in Table 5 for the results of the two different analy- 
sis methods and for various values for /x and n. It can be seen, that the Bayes 
method reproduces even the smallest lines properly, while the Maximum Like- 
hood method has some limitations there. 
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Table 5. Results from the analysis of the simulated spectra using a mean background 
of 4 and 0.5 counts, and different line intensities. The number d of cases, where the 
true number of counts in the line (given in the left column), is found in the calculated 
confidence area is given in the second or third column for the Bayes method, and in 
the fourth and fifth column for the maximum likelihood method. For details see text. 



Another test has been performed. We generated 1000 simulated spectra con- 
taining no line. Then the probability has been calculated with the Bayesian 
method that the spectrum does contain a line at a given energy. Table 6 presents 
the results: the first column contains the corresponding confidence limit c (pre- 
cisely the parameter defined earlier), the second column contains the expected 
number of spectra indicating existence of a line with a confidence limit above the 
value c and the third column contain the number of spectra with a confidence 
limit above the value c, found in the simulations. The result underlines that 
here is equivalent to the usual confidence level of classical statistics. 

We further investigated with the computer-generated spectra the dependence 
of the peak analysis on the width of the energy range of evaluation. Two examples 
are shown in Fig. 9. Here the contents of the simulated peak found with the 
Bayes method is shown as function of the analysis interval given in channels (one 
channel corresponds to 0.36 keV in our measured spectra). The line in the middle 
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C.L. 


Expected 


Found 


90.0% 


100 ± 31 


96 


95.0% 


50 ± 7 


42 


99.0% 


10 ± 3 


12 


99.9% 


1 ± 1 






Table 6. Number of spectra with a calculated confidence limit above a given value. 
For details see text. 



is the best-fit value of the method, the upper and lower lines correspond to the 
upper and lower 68.3% confidence limits. In the upper figure the true number of 
counts in the simulated line was 5 events, on a Poisson-distributed background 
of 0.5 events/channel, in the lower figure it was 20 events on a background of 
4 events/channel. It can be seen, that the analysis gives safely the correct number 
of counts, when choosing an analysis interval of not less than 40 channels. 

3.3 Analysis of the Full Data 

We first concentrate on the full spectra (see Fig. 4, and Fig. 2 in [1]), without 
any data manipulation (no background subtraction, no pulse shape analysis). 
For the evaluation, we consider the raw data of the detectors. 

The Bayesian peak finding procedure described in the last section leads to 
the result shown on the left hand sides of Figs. 11,12. For every energy E of 
the spectral range 2000 -2080 keV, we have determined the probability K e that 
there is a line at E. All the remainder of the spectrum was considered to be 
background in this search. 

The peak detection procedure yields lines at the positions of known [79] weak 
7-lincs from the decay of 2U Bi at 2010.7, 2016.7, 2021.8 and 2052.9 keV. The 
lines at 2010.7 and 2052.9 keV are observed at a confidence level of 3.7 and 2.6cr, 
respectively. The observed intensities are consistent with the expectations from 
the strong Bi lines in our spectrum, and the branching ratios given in [79] , within 
about the 2a experimental error (see Table 7 and [109]). The expectations here 
are calculated including summing effects, by Monte-Carlo simulation of our set- 
up. Only in this way the strong dependence of the relative intensities on the 
location of the impurities in the set-up can be properly taken into account (see 
Fig. 10). (A separate measurement with a g§ 6 -ffa source being in progress, will 
allow to study the intensities of the weak 214 Bi lines in the setup with high 
statistics). 

In addition, a line centered at 2039 keV shows up. This is compatible with 
the Q-value [77,78] of the double beta decay process. We emphasize, that at this 
energy no 7-line is expected according to the compilations in [79] (see discussion 
in section 2). Figs. 11,12 do not show indications for the lines from 56 Co at 
2034.7 keV and 2042 keV discussed earlier [76] (but sec also Fig. 14). We have 
at present no convincing identification of the lines around 2070 keV indicated by 
the peak identification procedure. 
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Fig. 9. Upper part: Analysis ol simulated spectrum with Gaussian peak of 5 events 
and FWHM of 9.4 channels on a Poisson-distributed background spectrum of 
0.5 events/channel, as function of interval of analysis. The middle line is the best value. 
Upper and lower lines correspond to the 68.3% confidence limits. Lower part: The same 
as above, but the peak contains 20 events, the background is 4 events/channel. One 
channel corresponds to 0.36 keV in our measured spectra. 



It may be important to note, that essentially the same lines as found by the 
peak scan procedure in Figs. 11,12,13, are found [107,109] when doing the same 
kind of analysis with the best existing natural Ge experiment of D. Caldwell et 
al. [110], which has a by a factor of ~4 better statistics of the background. This 
experiment does, however, not see the line at 2039 keV (see section 3.5). 

Bayesian peak detection (the same is true for Maximum Likelihood peak 
detection) of our data suggests a line at Qpp whether or not one includes detector 
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E(keV) 



Fig. 10. Demonstration of the dependence of the relative intensities of the weak Bi 
lines on the location of the impurity in the experimental setup. The right side shows 
two different locations of the impurity (1) relative to the detectors. The left side shows 
the corresponding different relative intensities to be expected. Upper part: Source in 
some distance from the detectors (as in a) right part). Lower part: Source in copper 
cap very near around the detectors (as in b) right part). The ordinate shows relative 
intensities (see [106]). 

Nr. 4 without muon shield (Figs. 11,12). The line is also suggested in Fig. 13 
after removal of multiple site events (MSE), see below. 

On the left-hand side of Figs. 11,12,13, the background intensity (l-?7) iden- 
tified by the Bayesian procedure is too high because the procedure averages the 
background over all the spectrum (including lines) except for the line it is try- 
ing to single out. Inclusion of the known lines into the fit naturally improves 
the background. As example, we show in Fig. 14 the spectrum of Fig. 4 (here 
in the original 0.36 keV binning) with a simultaneous fit of the range 2000 - 
2060 keV (assuming lines at 2010.78, 2016.70, 2021.60, 2052.94, 2034.76, 2039.0, 
2041.16keV). The probability for a line in this fit at 2039 keV is 86%. 
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Table 7. Bi is product of the U natural decay chain through /3~ decay of Pb 
and a decay of 218 At. It decays to 214 Po by (3~ decay. Shown in this Table are the 
measured intensities of 214 Bi lines in the spectrum shown in Fig.l of Ref. [1] in the 
energy window 2000 - 2060 keV, our calculation of the intensities expected on the basis 
of the branching ratios given in Table of Isotopes [79] , with and without simulation of 
the experimental setup, and the intensities expected by Aalseth et al. hep-ex/0202018, 
who do not simulate the setup and thus ignore summing of the 7 energies (see [4]). 
*) We have considered for comparison the 3 strongest 214 Bi lines, leaving out the line 
at 1120.287 keV (in the measured spectrum this line is partially overimposed on the 
1115.55 keV line of 65 Zn). The number of counts in each line have been calculated by a 
maximum-likelihood fit of the line with a gaussian curve plus a constant background. 
+) The simulation is performed assuming that the impurity is located in the copper 
part of the detector chamber (best agreement with the intensities of the strongest 
lines in the spectrum). The error of a possible misplacement is not included in the 
calculation. The number of simulated events is 10 s for each of our five detectors. 
**) This result is obtained normalizing the simulated spectrum to the experimental one 
using the 3 strong lines listed in column one. Comparison to the neighboring column 
on the right shows that the expected rates for the weak lines can change strongly if 
we take into account the simulation. The reason is that the line at 2010.7 keV can 
be produced by summing of the 1401.50 keV (1.55%) and 609.31 keV (44.8%) lines, 
the one at 2016.7 keV by summing of the 1407.98 (2.8%) and 609.31 (44.8%) lines; 
the other lines at 2021.8 keV and 2052.94 keV do suffer only very weakly from the 
summing effect because of the different decay schemes. 

++) This result is obtained using the number of counts for the three strong lines 
observed in the experimental spectrum and the branching ratios from [79] , but without 
simulation. For each of the strong lines the expected number of counts for the weak 
lines is calculated and then an average of the 3 expectations is taken. 
***) Without simulation of the experimental setup. The numbers given here are close 
to those in the neighboring left column, when taking into account that Aalseth et al. 
refer to a spectrum which has only ~11% of the statistics of the spectrum shown in 
Fig.l of Ref. [1]. 
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Fig. 11. Scan for lines in the full spectrum taken from 1990-2000 with detectors Nr. 
1,2,3,4,5, (Fig. 4), with the Bayesian method of sec. 3.1. The ordinate is the probability 
Ke that a line exists as defined in the text. Left: Energy range 2000 - 2080 keV. 
Right: Energy range of analysis ±5<r around Qpp. 
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Fig. 12. Left: Probability Ke that a line exists at a given energy in the range of 2000- 
2080 keV derived via Bayesian inference from the spectrum taken with detectors Nr. 
1,2,3,5 over the period August 1990 to May 2000, 46.502 kgy (see Fig.2 from [1].) Right: 
Result of a Bayesian scan for lines as in the left part of this figure, but in the energy 
range of analysis ±5cr around Q/3/3. 
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Fig. 13. Scan for lines in the single site event spectrum taken from 1995-2000 with 
detectors Nr. 2,3,5, (Fig. 5), with the Bayesian method (as in Figs. 11,12). Left: Energy 
range 2000 -2080 keV. Right: Energy range of analysis ±4.4u around Q/3/3. 



Finally, on the right-hand side of Figs. 11,12 (and also Fig. 13) the peak detection 
procedure is carried out within an energy interval that seems to not contain 
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Fig. 14. Simultaneous analysis of the spectrum measured with the Ge detectors Nr. 
1,2,3,4,5 over the period August 1990 - May 2000 (54.9813 kgy) (same as in Fig. 4, 
but here shown in the 0.36 keV original binning) in the energy range 2000 - 2060 keV, 
with the Maximum Likelihood Method. The probability for a line at 2039.0 keV found 
is this way is 86%. 



(according to the left-hand side) lines other than the one at Q/3/3. This interval 
is broad enough (about ± 5 standard deviations of the Gaussian line, i.e. as 
typically used in search for resonances in high-energy physics) for a meaningful 
analysis (see Fig. 9 in section 3.2). We find, with the Bayesian method, the 
probability K# = 96.5% that there is a line at 2039.0 keV is the spectrum shown 
in Fig. 4. This is a confidence level of 2.1 a in the usual language. The number of 
events is found to be 0.8 to 32.9 (7.6 to 25.2) with 95% (68%) c.L, with best value 
of 16.2 events. For the spectrum shown in Fig. 2 in [1], we find a probability for 
a line at 2039.0 keV of 97.4% (2.2 a). In this case the number of events is found 
to be 1.2 to 29.4 with 95% c.L. It is 7.3 to 22.6events with 68.3% c.L. The most 
probable number of events (best value) is 14.8 events. These values are stable 
against small variations of the interval of analysis, as expected from Fig. 9 in 
section 3.2. For example, changing the lower and upper limits of the interval of 
analysis between 2030 and 2032 and 2046 and 2050 yields consistently values of 
K E between 95.3 and 98.5% (average 97.2%) for the spectrum of Fig.2 of [1]. 

We also applied the Feldman-Cousins method recommended by the Particle 
Data Group [88,89]. This method (which does not use the information that the 
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line is Gaussian) finds a line at 2039 keV on a confidence level of 3.1 a (99.8% 
c.L). 

3.4 Analysis of Single Site Events Data 

From the analysis of the single site events (Fig. 5), we find after 28.053 kg y of 
measurement 9 SSE events in the region 2034.1 - 2044.9 keV (± 3tr around Q^) 
(Fig. 15). Analysis of the single site event spectrum (Fig. 5), as described in 




Fig. 15. Events classified to be single site events (SSE) by all three methods of PSA in 
the range 2034.1 - 2044.9 keV, in the measurement period 10.1995 - 05.2000, 28.053 kgy. 

section 3.1, shows again evidence for a line at the energy of Qpp (Fig. 13). Ana- 
lyzing the range of 2032 - 2046 keV, we find the probability of 96.8% that there 
is a line at 2039.0 keV. We thus see a signal of single site events, as expected for 
neutrinoless double beta decay, precisely at the value obtained in the preci- 
sion experiment of [77]. The analysis of the line at 2039.0 keV before correction 
for the efficiency yields 4.6 events (best value) or (0.3 - 8.0) events within 95% 
c.L ((2.1 - 6.8) events within 68.3% c.L). Corrected for the efficiency to identify 
an SSE signal by successive application of all three PSA methods, which is 0.55 
± 0.10, we obtain a Ov/3/3 signal with 92.4% c.L. The signal is (3.6 - 12.5) events 
with 68.3% c.L (best value 8.3 events). Thus, with proper normalization con- 
cerning the running times (kgy) of the full and the SSE spectra, we see that 
almost the full signal remains after the single site cut (best value), while the 
21A Bi lines (best values) are considerably reduced. The reduction is comparable 
to the reduction of the 2103 keV and 2614 keV 228 Th lines (known to be multiple 
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site or mainly multiple site), relative to the 1592 keV 22S Th line (known to be 
single site), see Fig. 16. The same reduction is also found for the strong 214 Bi 
lines (e.g. at 609.6 and 1763.9keV (Fig. 16)). 




Fig. 16. Relative suppression ratios: Remaining intensity after pulse shape analysis 
compared to the intensity in the full spectrum. Right: Result of a calibration measure- 
ment with a Th source - ratio of the intensities of the 1592 keV line (double escape 
peak, known to be 100% SSE), set to 1. The intensities of the 2203 keV line (single 
escape peak, known to be 100% MSE) as well as of the 2614 keV line (full energy peak, 
known to be dominantly MSE) are strongly reduced (error bars are ±lcr). The same 
order of reduction is found for the strong Bi lines occurring in our spectrum (see Fig. 3) 
- shown in this figure are the lines at 609.4 and 1763.9 keV. Left: The lines in the range 
of weak statistics around the line at 2039 keV (shown are ratios of best fit values). The 
Bi lines are reduced compared to the line at 2039 keV (set to 1), as to the 1592 keV 
SSE Th line. 



The Feldman- Cousins method gives a signal at 2039.0 keV of 2.8 a (99.4%). 
The possibility, that the single site signal at 2039.0 keV is the double escape 
line corresponding to a (much more intense!) full energy peak of a 7-line at 
2039+ 1022=3061 keV is excluded from the high-energy part of our spectrum 
(see Fig. 17). 

3.5 Comparison with Earlier Results 

We applied the same methods of peak search as used in sections 3.3, 3.4, to the 
spectrum measured in the Ge experiment by Caldwell et al. [110] more than a 
decade ago. These authors had the most sensitive experiment using natural Ge 
detectors (7.8% abundance of 76 Ge). With their background being a factor of 
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Fig. 17. The measured sum spectrum up to ~ 8MeV (upper part) (22.36 kgy). No 
lines can be identified in the range 3-8 MeV, except a peak at 5.24±0.03 MeV, which 

Lower part: the range 3000 - 



we identify as a decay from 210 Po to 206 Pb (see [76]) 



3100 keV, where the full energy peak should occur at 3061 keV, in case the 2039 keV 
signal would be the double escape peak of a 7-line. 



9 higher than in the present experiment, and their measuring time of 22.6 kgy, 
they have a statistics for the background larger by a factor of almost 4 in their 
(very similar) experiment. This allows helpful conclusions about the nature of 
the background. 

The peak scanning finds (Fig. 18) indications for peaks essentially at the same 
energies as in Figs. 11,12,13. This shows that these peaks are not fluctuations. 
In particular it sees the 2010.78 and 2052.94 keV 214 Bi lines with 3.6 and 2.8 
g c.L, respectively. It finds, however, no line at Qpp (see also Fig. 19). This is 
consistent with the expectation from the rate found from the HEIDELBERG- 
MOSCOW experiment. About 17 observed events in the latter correspond to 0.6 
expected events in the Caldwell experiment, because of the use of non-enriched 
material and the shorter measuring time. 

Another Ge experiment (IGEX) using 9 kg of enriched 76 Ge, but collecting 
since beginning of the experiment in the early nineties till shutdown in end of 
1999 only 8.8 kgy of statistics [114], because of this low statistics also naturally 
cannot see any signal at 2039 keV. 
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Fig. 18. Peak scanning of the spectrum measured by Caldwell et al. [110] with the 
Maximum Likelihood method (upper part), and with the Bayesian method (lower part) 
(as in Fig. 11,12,13) (see [109]). 



3.6 Some Comments on the Bayesian, x 2 an d 
Maximum Likelihood Methods 

We probed the sensitivity of peak identification for the three methods: Bayesian, 
X 2 , Maximum-Likelihood, for the latter two using codes from [92]. 

The disadvantage of the latter two methods is, that at low counting rates 
observation of lines with negative counting rate is possible. This is excluded in 
the Bayesian method. We find that the Bayesian method tends to systematically 
give too conservative confidence limits (see Table 5 in section 3.2). We shall 
discuss technical details of the three methods in a separate paper [109]. 



4 Half-life of the Neutrinoless Mode and 
Effective Neutrino Mass 



We emphasize that we find in all analyses of our spectra a line at the value 
of Q/3/3- We have shown that to our present knowledge the signal at does 
not originate from a background 7-line. On this basis we translate the observed 
number of events into half-lives for the neutrinoless double beta decay. We give 
in Table 8 conservatively the values obtained with the Bayesian method and not 
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2050 2060 
energy(KeV) 



Fig. 19. Analysis of the spectrum measured by D. Caldwell et al. [110], with the 
Maximum Likelihood Method, in the energy range 2000 - 2060 keV (as in Fig. 14) 
assuming lines at 2010.78, 2016.70, 2021.60, 2052.94, 2039.0 keV. In contrast to Fig. 14 
no indication for a signal at 2039 kev is observed (see text) . 



those obtained with the Feldman-Cousins method. Also given in Table 8 are the 
effective neutrino masses (m) deduced using the matrix elements of [33,32]. 

We derive from the data taken with 46.502 kg y the half-life T°J 2 = (0.8 - 
18.3) x 10 25 y (95% c.l.). The analysis of the other data sets, shown in Table 8 
confirm this result. Of particular importance is that we see the 0v(3/3 signal in 
the single site spectrum. 

The result obtained is consistent with the limits given earlier by the HEI- 
DELBERG-MOSCOW experiment [19]. It is also consistent with all other double 
beta experiments - which still reach less sensitivity (see Figs. 21,20). A second 
Ge-experiment [114], which has stopped operation in 1999 after reaching a sig- 
nificance of 9kgy, yields (if one believes their method of 'visual inspection' in 
their data analysis) in a conservative analysis a limit of T°^ 2 > 0.55 x 10 25 y 
(90% c.l.). The 128 Te geochemical experiment [111] yields {m v ) < 1.1 eV (68 % 
c.l.), the 130 Te cryogenic experiment yields [117] (m„) < 1.8 eV and the CdW04 
experiment [112] (m„) < 2.6 eV, all derived with the matrix elements of [33] to 
make the results comparable to the present value. 

Concluding we obtain, on the above basis, with more than 95% probability, 
first evidence for the neutrinoless double beta decay mode. 

As a consequence, at this confidence level, lepton number is not conserved. 
Further our result implies that the neutrino is a Majorana particle. Both of these 
conclusions are independent of any discussion of nuclear matrix elements. 

The matrix element enters when we derive a value for the effective neutrino 
mass. If using the nuclear matrix element from [33,32], we conclude from the 
various analyses given above the effective mass (m) to be (m) = (0.11 - 0.56) eV 
(95% c.l.), with best value of 0.39 eV. Allowing conservatively for an uncertainty 
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Fig. 21. Present evidence for 0v/3f3 decay from data of the HEIDELBERG-MOSCOW 
experiment and the potential of present and future (3/3 experiments. Vertical axis - 
effective neutrino mass limits (only HEIDELBERG-MOSCOW gives a value) in eV, 
horizontal - isotopes used in the various experiments. 



of the nuclear matrix elements of ± 50% (for detailed discussions of the status 
of nuclear matrix elements we refer to [25,26,34,37,6,28,42,40,41]) this range 
may widen to (m) = (0.05 - 0.84) eV (95% c.l.). In Table 9 we demonstrate the 
situation of nuclear matrix elements by showing the neutrino masses deduced 
from different calculations. It should be noted that the value obtained in Large 
Scale Shell Model Calculations [35] is understood to be too large by almost a 
factor of 2 because of the two small configuration space, (see, e.g. [37]), and that 
the second highest value given (from [36]), has now been reduced to 0.53 eV [42]. 
The recent studies by [40,41] yield an effective mass of (0.44 - 0.52) eV. We see 
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Significan- 
ce [kgy] 


Detectors 




(m) eV 


Conf. 
level 


54.9813 


1,2,3,4,5 


(0.80 - 35.07) x 10 25 
(1.04 - 3.46) x 10 25 
1.61 x 10 25 


(0.08 - 0.54) 
(0.26 - 0.47) 
0.38 


95% c.l. 
68% c.l. 
Best Value 


46.502 


1,2,3,5 


(0.75 - 18.33) x 10^ 
(0.98 - 3.05) x 10 25 
1.50 x 10 25 


(0.11 - 0.56) 
(0.28 - 0.49) 
0.39 


95% c.l. 
68% c.l. 
Best Value 


28.053 


2,3,5 SSE 


(0.88 - 22.38) x 10 2b 
(1.07 - 3.69) x 10 25 
1.61 x 10 25 


(0.10 - 0.51) 
(0.25 - 0.47) 
0.38 


90% c.l. 
68% c.l. 
Best Value 



Table 8. Half-life for the neutrinoless decay mode and deduced effective neutrino mass 
from the HEIDELBERG-MOSCOW experiment. 



M°" from 


[32,33] 


[34] 


[31] 


[35] 


[36] 


[38] 


[42] 


[40,41] 


(m) cV 95% 


0.39 


0.37 


0.34 


1.06 


0.87 


0.60 


0.53 


0.44 - 0.52 



Table 9. The effect of nuclear matrix elements on the deduced effective neutrino 
masses. Shown are the neutrino masses deduced from the best value of T*jy 2 = 1.5 x 
10 25 y determined in this work with 97% c.l. when using matrix elements from various 
calculations and a phase factor of F°" — 6.31 x 10~ 15 y _1 . 



that the early calculations [33] done in 1989 agree within less than 25% with the 
most recent values. 

In the above conclusion for the effective neutrino mass, it is assumed that 
contributions to 0^/3/3 decay from processes other than the exchange of a Majo- 
rana neutrino (see, e.g. [6,94] and section 1) are negligible. It has been discussed, 
however, recently [55] that the possibility that Of/3/3 decay is caused by R-parity 
violation, may experimentally not be excluded, although this would require mak- 
ing R-parity violating couplings generation-dependent. 

5 Conclusions and Outlook 

With the value deduced for the effective neutrino mass, the HEIDELBERG- 
MOSCOW experiment excludes several of the neutrino mass scenarios allowed 
from present neutrino oscillation experiments (see Fig. 22) - allowing mainly 
only for degenerate mass scenarios, and an inverse hierarchy iv and Av- scenario 
(the former of these being, however, strongly disfavored by a recent analysis 
of SN1987A [47]). For details we refer to [3]. In particular, hierarchical mass 
schemes are excluded at the present level of accuracy. 
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Fig. 22. The impact of the evidence obtained for neutrinoless double beta decay in this 
paper (best value of the effective neutrino mass (m) = 0.39 eV, 95% confidence range 
(0.05 - 0.84) eV - allowing already for an uncertainty of the nuclear matrix element 
of a factor of ± 50%) on possible neutrino mass schemes. The bars denote allowed 
ranges of (m) in different neutrino mass scenarios, still allowed by neutrino oscillation 
experiments. Hierarchical models are excluded by the new 0vf3j3 decay result. Also 
shown are the expected sensitivities for the future potential double beta experiments 
CUORE, MOON, EXO and the 1 ton and 10 ton project of GENIUS [96] (from [3]). 



According to [51] a global analysis of all solar neutrino data including the 
recent SNO neutral-current rate selects the Large Mixing Angle (LMA) at the 
90% c.L, however, the LOW solution is also viable, with 0.89 goodness of fit. 

Assuming the degenerate scenario to be realized in nature we fix - according 
to the formulae derived in [5] - the common mass eigenvalue of the degenerate 
neutrinos tom= (0.05 - 3.4) eV. Part of the upper range is already excluded by 
tritium experiments, which give a limit of m < 2.2 eV, or 2.8 eV (95% c.l.) [95]. 
The full range can only partly (down to ~ 0.5 eV) be checked by future tritium 
decay experiments, but could be checked by some future (3(3 experiments (see, 
e.g. [6,96]). The deduced 95% interval for the sum of the degenerate neutrino 
masses is consistent with the range for Q v deduced from recent cosmic microwave 
background measurements and large scale structure (redshift) surveys, which 
still allow for a J2i m i — 4.4 eV [104,103]. The range of (m) fixed in this work 
is, already now, in the range to be explored by the satellite experiments MAP 
and PLANCK [102,105] (see Fig. 23). It lies in a range of interest for Z-burst 
models recently discussed as explanation for super-high energy cosmic ray events 
beyond the GKZ-cutoff [99,100,66]. Finally, the deduced best value for the mass 
is consistent with expectations from experimental /i — > branching limits 
in models assuming the generating mechanism for the neutrino mass to be also 
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m,(eV) HEIDELBERG-MOSCOW Positive EVIDENCE 




Fig. 23. Double beta decay observable (m) and oscillation parameters: The case for 
degenerate neutrinos. Plotted on the axes are the overall scale of neutrino masses mo 
and mixing tan 2 8 12 . Also shown is a cosmological bound deduced from a fit of CMB 
and large scale structure [102] and the expected sensitivity of the satellite experiments 
MAP and PLANCK. The present limit from tritium j3 decay of 2.2 eV [95] would lie 
near the top of the figure. The range of (m) fixed by the HEIDELBERG-MOSCOW 
experiment is, in the case of small solar neutrino mixing, already in the range to be 
explored by MAP and PLANCK. 

responsible for the recent indication for an anomalous magnetic moment of the 
muon [97]. A recent model with underlying A4 symmetry for the neutrino mixing 
matrix (and the quark mixing matrix) also leads to degenerate neutrino masses 
consistent with the present experiment [98] . This model succeeds to consistently 
describe the large (small) mixing in the neutrino (quark) sector. 

The neutrino mass deduced leads to 0.002 < Sl^h 2 < 0.1, and thus may allow 
neutrinos to still play an important role as hot dark matter in the Universe (see 
also [54]). 

With the HEIDELBERG-MOSCOW experiment, the era of the small smart 
experiments is over. New approaches and considerably enlarged experiments (as 
discussed, e.g. in [6,29,96,101]) will be required in future to fix the neutrino mass 
with higher accuracy. 

Since it was realized in the HEIDELBERG-MOSCOW experiment, that the 
remaining small background is coming from the material close to the detector 
(holder, copper cap, ...), elimination of any material close to the detector will 
be decisive. Experiments which do not take this into account, like, e.g. CUORE 
[117], and MAJORANA [115], will allow only rather limited steps in sensitivity. 

Another crucial point is - see eq. (6) - the energy resolution, which can be 
optimized only in experiments using Germanium detectors or bolometers. It will 
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be difficult to probe evidence for this rare decay mode in experiments, which 
have to work - as result of their limited resolution - with energy windows around 
Q/3/3 of up to several hundreds of keV, such as NEMO III [118], EXO [116], 
CAMEO [113]. 

For example - according to eq. (6) - to compensate for the projected energy 
resolution of only 130 keV of EXO [116], the mass of the EXO experiment has 
to be increased to almost half a ton of enriched material, to reach the sensitivity 
of the HEIDELBERG-MOSCOW experiment. Only after that one can think 
about improving the sensitivity, by improving the background. Correspondingly, 
according to eq. (6), a potential future 100 kg 82 Se NEMO experiment would be 
because of its low efficiency, equivalent only to a 10 kg experiment (not talking 
about the energy resolution) . 

In the first proposal for a third generation double beta experiment, the GE- 
NIUS proposal [21,96], the idea is to use 'naked' Germanium detectors in a huge 
tank of liquid nitrogen. It seems to be at present the only proposal, which can 
fulfill both requirements mentioned above. The potential of GENIUS is together 
with that of some later proposals indicated in Fig. 22. GENIUS would - with only 
100 kg of enriched 76 Ge - increase the confidence level of the present 0v[3f3 signal 
to 5cr within one year of measurement. A GENIUS Test Facility is at present 
under construction in the GRAN SASSO Underground Laboratory [119]. 
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